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Abstract. The derivation of kinetic equation for one-particle distribution function F{p,r,t) from 
BBGKY chain leads to the collision integral expressed in terms of the two-particle correlation 
function. The latter is in turn expressed by F using Boltzmann's Stosszahlansatz that implies 
the neglect of initial one-time correlation function g2{pl,p2,r,t) which is non-diagonal in the 
momentum space. In non-equilibrium case, however, it has been established that pair collisions 
generate the non-diagonal two-particle correlations that give the essential contribution to the current 
fluctuations of hot electrons. These correlations have been shown giving also a contribution to the 
collision integrals, i.e., to kinetic properties of nonequilibrium gas The expression for electron 
energy loss rate P via phonons is re-derived in detail from this point of view. The order of value 
of the new term in P in the cases of acoustic- and optic al-phonon scattering is evaluated in the 
electron temperature approximation using the experimental study of hot-electron energy relaxation 
time in « — InSb at helium temperatures. 
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INTRODUCTION 

It was established that pair collisions generate in non-equilibrium case non-zero two- 
particle correlations that are non-diagonal in momentum space and give the essential 
contribution to the current fluctuations of hot electrons These correlations have 

been shown giving also a contribution to the collision integrals, i.e., to kinetic proper- 
ties of non-equilibrium gas IHH]. The expression for electron energy loss rate P via 
phonons is re-derived from this point of view. The formula for P in terms of electron 
density-density correlation function {pkPk) a was obtained by Kogan [5] using Fermi 
Golden Rule. However, in order to come to the well-known result for P the fluctuation- 
dissipation theorem was implicated to express {pkp^) ^ i'^ terms of the dielectric func- 
tion of the electron gas. On this way the non-equilibrium correction has been lost. 
Here the contribution of the non-diagonal part of the non-equilibrium pair correlator 
to phonon-electron collision integral and to P is obtained and explicitly calculated in 
the electron temperature approximation. The contribution obtained does not contain the 
extra powers of the small parameter of kinetic theory unlike the virial expansions in 
equilibrium thermodynamics. The order of value of the new term in P is evaluated us- 
ing the experimental study of hot-electron energy relaxation time in n — InSb at helium 
temperatures |6]. 



BASIC DESIGNATION 



Hamiltonian 

H = He+Hp + Vep 

Electron part 

He = L^Qk^k + Vee, 
k 

where V^e stands for the Coulomb electron-electron interaction 

Phonon part 

Hp = l^hCDtb^bt 
f 

and electron-phonon interaction 

Vep^^^i:{cfP+bt + c}Ptb^), 

where Vq is the normalizing volume. 

Fourier components of the electron density operator 

^ k+2 5 

Phonon distribution function 

A^f(0 = ((^f+^f),> 

Electron distribution function 

Fk(0 = (Kak),> 

The angle brackets denote the averaging with the density matrix corresponding to the 

noninteracting electrons and phonons at initial instant of time. The subscript t indicates 
the Heisenberg representation for electron and phonon operators. 



KINETIC EQUATION FOR PHONONS 

We obtain the kinetic equation for slow variables Nf{t) , smoothing Heisenberg equa- 
tions of motion over mesoscopic time interval. With accuracy to the second-order terms 
ine — p interaction we have 

/ t+At tl \ 

^ = (^^^(0 / dtt I dtn {btb,)^\] S{t) 

The superscript denotes operators in the interaction representation with respect to the 
electron-phonon interaction Vgp . 

Calculations of the internal commutator lead to the desired result involving the elec- 
tron density-density correlation function: 

where the condition (Of A? >> 1 was used. The function 

oo 

{ptPt')co= J dT{pt{t+T)p + {t))QXp{i(OT) 



is the spectral density of the correlation function of the electron density fluctuations. 



Spectral density of the high-frequency electron fluctuations in the 
electron temperature approximations 

The spectral density must be calculated without Vep taken into account. From equa- 
tions of motion for the electron operators we obtain in the random phase approximation: 

where e((D, f) is the relative dielectric function of the electron gas embedded into the di- 
electric medium and (Pf^Pf)^ is the spectral density of the density-density correlation 
function for non-interacting electrons. The latter is connected with one-time correla- 
tion function ^Pf(k)Pj^(ki))^ by equations of motion for free electrons. The one-time 
correlation function is defined as 

(Pf(k)p+(ki))^ = (4_f«k+|<, + f%-i/^ 

The one-time correlation function is independent of time (in steady state) under follow- 
ing condition on the electron energies = hk^ /2m : 

^k-i - ^k+l + ^k,+f - ^k,-| = (ki - k) = 0. 

Because f 7^ then it is necessary to be ki = k . Therefore, this condition provides 
the one-time correlator to be the integral of motion during the time interval A? like the 
distribution functions A'f and F^. 



Kinetic equation with account for the contribution of the 
non-equilibrium pair correlations 

The substitution of the explicit expression Eq.(0)for correlators to the kinetic equation 
Eq.O gives rise to 







dt 


Von^\£{(Oi,t)\^ 




2|C/|' 



yo^1e(«f,f)| k 



L k-f k+f )^ ^ (^k+f - ^k-l - ^Of) (A^f + 1) - 

(3) 

E k+i<+f k-i ) ^ (^k+i - ^k- 1 - «f) ^f- 



Taking into account the commutation rules for a^a^^ it is easy to obtain 



According to the results from [1] 



% «ki «k ) = -Pk^k, + (k, ki ; 



where ^ (k, ki) is a non-equilibrium correction to the pair correlation function. 

Finally, the kinetic equation for phonons with non-equilibrium corrections of electron 
pair correlators takes the form: 



dNf 



2 C 



f\ 



dt yo^^|e(%,f)|^ 

xl[^k+| (l-^k4)-'^>(k + ^,k-0](A^f+l)5(^2,^,-^2,_,-cOf) 

(l -i^+i) -<l>(k-f/2,k + f/2)] Ntd - (Of) 



(4) 



Certainly, it needs to keep in mind the Eq.© describes the part of full time dependence 
of the phonon occupation numbers that resulting from electron-phonon interaction only. 

It should be stressed that for ordinary classical equilibrium electron gas there are no 
pair correlations between different states of momentum space. This fact can be easily 
seen from the form of the Gibbs distribution. In conclusion of these considerations 
it needs to underline that the contribution of non-equilibrium correlations to kinetics 
should be expected in the cases when the particles considered are scattered by non- 
equilibrium many-particle system. 



ELECTRON ENERGY LOSS RATE VIA PHONONS 

The energy loss rate of electrons due to phonons per unit volume can be obtained using 
the relation 




It is seen from the kinetic equation for phonons Eq.® that the loss rate can be broken 
up into two parts: 

P = Pi+P2, 

where the first term on the right hand side is related to the one-particle electron distri- 
bution function and it is described by the known expressions. The second one is due to 
the correlation contribution determined by the function .We assume during the calcu- 
lations of P the electron-electron collisions are sufficiently frequent in order to provide 



the use of the electron temperature approximation for hot electrons. As a result of these 
assumptions the following formulae may be derived 



f {Voh)^\£{cOf,{)\ 



X 



k 

and 



2\Cf 



2 



P2{T, To) = - £ „2 ^^E (k + f , k) 5 (Qk+f - - G>f ) (6) 

Here Nf{T) = [exp{h(Of/T) — 1]^^ is the equilibrium Bose-Einstein distribution with 
electron temperature T and A^f is the true non-equilibrium phonon distribution function. 
Further, we neglect possible phonon heating and consider that the phonons remain in 
thermal equilibrium at temperature Tq, so Nf = (Tq). 

The function ^ (k, ki) was found in |2] under accepted assumptions for the case of 
nondegenerate electron gas and can be presented in the form: 

dF^ dF^ 

0(k,ki) = e(r,ro)^^ (7) 



where 

Q{T,To) = ^ 

nce 



P{TJq)/{T-Tq) 
dP/dT 



-1 



Here Nc{T) and nCe = 3n/2 are the effective density of states and specific heat capaci- 
tance of electron gas with concentration n, respectively. The expression for in ^\ was 
obtained without accounting for the contribution of itself in P .Therefore, the quantity 
P\ was only implied as the total energy loss rate. Now evidently we have the differential 
equation for determination of P2{T,Tq) or full P since this quantity is defined by the 
right hand side of the equality Eq.© including P2{T,Tq) itself. The situation is quite 
resembling the self-consistent field approach in plasma theory. It needs to stress that the 
substitution of the total loss rate P instead of Pi in the expression for (k, ki ) is justified 
because all symmetry properties of the electron-electron and electron-phonon integrals 
of collisions used in the derivation are not destroyed by introducing (k, ki ) in collision 
integrals. 

It is of interest to estimate the interrelation between P\ and P2- Instead carrying out 
the integrations in Eqs. (|5l)-® and following solution of a nonlinear differential equation 
for P (or P2) we can use the experimental data from Ref. [6] where hot-electron energy 
relaxation time at helium temperatures in n — InSb was investigated and all necessary 
information on n, P, dP/dT, and T can be extracted from. It was shown there that at 
low heating the electron energy loss is due to acoustic phonon scattering and next it is 
changed by longitudinal optical (LO) phonon scattering. So we can evaluate P2/P1 for 
both types of the electron-phonon scattering. 



It is necessary to carry out preliminary the integrations over k in Eqs.©-® for all 
essential parameters be explicitly introduced. Taking into account nondegeneracy of 
electron gas the expressions for Pi and P2 can be transform as following 



P,{TJq) = {T-Tq) 



2ms 



2 Cf 



{27t)^h^ y hVo\£{(Of,{ 
9Tm^s \P{TJq)/{T-To) . 



|2^f/2 



NciT) 



exp 



dP/dT 



2mT 



2 Cf 



t hVo\e{(Of,{)\ 



(8) 



yF,%iT/2)[\ + 0{T)] 
(9) 



for acoustic phonon scattering. Here s is the sound velocity and it has been accounted for 
the near-elastic scattering of electrons by acoustic phonons. While the sign of expression 
I© for P2 leaves indeterminate since implicit 0{T) terms may redefine it, the order of 
magnitude of P2/P1 can be found. The expression Eq.® is just the Kogan's formula 
from ijsl] corrected for the screening of electron-phonon interaction. 
In the case of LO phonon scattering we obtain 



^1(7^, 7b) 



m TficOto ^ 



X 



' NciT) 
2\Cf\^ 



2 



f nvo\£io}f,f)\ f 



■exp(- 



2mT 



LO 



2/ 



) 



(10) 



P{TJq)/{T-Tq) 



2{2KYh^Nc{T) 



dP/dT 



exp 



hG)L0 \ 

T j 



exp(- 



2h' 



2mT 



LO 



2/ 



[\ + 0{T)] 



(11) 



where fj^g = 2ma)Lo/^ ■ 

The sums over f in corresponding expressions are of the same order of magnitude. 
Therefore, we have to evaluate following ratios 



Pi 
Pi 



T-To 

for the acoustic phonon scattering and 



P2 
Pi 



ST 



P{TJq)/{T-Tq) 
dP/dT 



P{T,To)/{T-To) 
dP/dT 



1 



(12) 



(13) 



for the low-temperature LO phonon scattering (h (Ow /T » 1 ) . Turning to experimental 
data 10] and using relation between dP/dT and measured electron energy relaxation 



time Tp = ncel {dP/dT) we obtain for the acoustic scattering P2/P1 ~ 1.3 at T ~ 1 .%K, 
P^6- IQ-^W/cm^, T/>~4-10-^5andP2M ^0.6atr~ 15.8if,/'^3.8-10-2w/cm3, 
Tp ~ 0.8 ■ 10^^ s for LO phonon scattering. All data are taken at Tq = 4.2K. These 
estimations may be considered as an indication that the pair correlations contribute to 
hot electron energy loss via phonons comparable with common accepted expressions. 



Concluding remarks 

It is of importance to note in conclusion that the expression for (k, ki ) found in 
the limit of high frequency of the electron-electron collisions does contain coupling 
constants of neither electron-electron nor electron-phonon interactions. Therefore, the 
corresponding corrections to the integral of phonon-electron collisions do not contain 
additional powers of the small parameters of the kinetic theory. The interesting but still 
open question: are there analogous contributions to the hydrodynamic equations? 
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